Lengthscale dependence of dynamic four-point susceptibilities in glass formers 
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Dynamical four-point susceptibilities measure the extent of spatial correlations in the dynamics 
of glass forming systems. We show how these susceptibilities depend on the length scales that 
necessarily form part of their definition. The behaviour of these susceptibilities is estimated by 
means of an analysis in terms of renewal processes within the context of dynamic facilitation. The 
analytic results are confirmed by numerical simulations of an atomistic model glass-former, and of 
two kinetically constrained models. Hence we argue that the scenario predicted by the dynamic 
facilitation approach is generic. 



I. INTRODUCTION 



The length scales governing dynamical heterogeneity 

in glass-forming liquids mail! 

are often described 
in terms of the susceptibility associated with fluctuations 
in the self intermediate scattering function [1, H, 0, H, Q : 



(i) 



Here the indices j and I run over the N particles in the 
system, the position of the jth particle at time t is fj(i); 



SF 3 {k,t) 



o ifc-[*»(t)-#j(0)] 



Jk-[r 3 (t)-r 3 {Q)] 



(2) 



and k — \k\. Under supercooled conditions, this four- 
point correlation function typically grows in time to- 
wards a peak, before decreasing at large times. This 
non-monotonic behaviour is a consequence of the tran- 
sient nature of dynamic heterogeneity. 

It was suggested by Toninelli el al. that the time 
dependence of \i{k,t) can be used to distinguish be- 
tween different theoretical scenarios for the glass tran- 
sition. The dependence of Xi(^,t) on the wave vector 
k was considered for a glass-formin g sy stem in [lOj, and 



for a sheared granular material in [11] . In both cases, 
significant dependence on wave vector was found. In this 
article, we investigate this wave vector dependence, in 
particular the way that Xi{k,t) grows towards its peak. 
We present data for an atomistic system, and for two 
kinetically constrained models [T3 |. We find non-trivial 
wave vector dependence in all three cases. We explain 
this generic behaviour analytically using a treatment that 
we used earlier to describe dynamic decoupling in glass 
formers (HQ El]. 

Our analysis shows that the non-trivial behaviour of 
four-point correlators comes from two sources. The first 
contribution arises because particles that have not moved 
are clustered in space; the second comes from correla- 
tions between particle displacements. The relative sizes 
of these contributions depend on the wave vector k. The 
first dominates when k is large, and the second dominates 
when k is small. The crossover between these two regimes 
corresponds to the crossover between non-Fickian and 
Fickian regimes observed in two-point functions [Tol l . 
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FIG. 1: (Color online) The distinct part of xi(k,t) as a func- 
tion of t for two wave vectors in three model systems: the 
supercooled three dimensional WCA mixture (a); a kineti- 
cally constrained triangular lattice gas (b) ; and probes in the 
one dimensional Fredrickson- Andersen model (c). The sym- 
bol xf U (k,t) denotes (\SFi(k, t)\ 2 }. Full details are presented 
in the relevant sections below. Dashed and dotted lines show 
that in all cases the exponents of power law fits increase with 
decreasing wave vector. The largest wave vectors considered 
are the principal wave vectors for each system (for the WCA 
case this is the peak location in the equilibrium structure fac- 
tor, g ). 



Figure [T] illustrates the behaviour that we consider for 
three different model systems. These systems, described 
in detail in sections IIIII and IIV1 represent three levels 
of coarse-graining in the glassy system. The most de- 
tailed is a fluid mixture of classical particles in contin- 
uous three dimensional space, interac ting with Weeks- 
Chandler- Andersen (WCA) potentials [laGjj]- The sec- 
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ond is the so-called (2)-TLG, due to Jackie and Kronig 
18]. It is a kinetically constrained lattice gas in which the 
dynamics are highly co-operative and relaxation times 
increase very quickly with increasing density. The fi- 
nal system is the one dimensional one spin facilitated 
Fredrickson-Andersen (FA) model to which probe 
particles have been added, following Ref . [l3| . The probe 
particles do not interact with one another, but they prop- 
agate through an environment that is dynamically het- 
erogeneous. In this model, all of the atomistic interac- 
tions have been removed, leaving only an idealised dy- 
namically heterogeneous system. 

The similarity of Xi(k,t) between these three different 
model systems is striking. The four-point susceptibility 
has the usual single peak, and the increasing part of each 
curve can be fitted by a power law. In all three cases, the 
exponent with which Xi(k,t) increases depends strongly 
on k. Further, all cases exhibit a shift of the peak of 
Xi{k,t) to later times as k decreases. 

II. THEORY OF FOUR-POINT FUNCTIONS IN 
HETEROGENEOUS SYSTEMS 

A. The dynamical facilitation approach 

Our starting point is to follow [TEj and write 

Fj(k,t) ee exp[ifc • Afj(t)} 

= Pj (t) + [1 - Pj (*)] exp [ik ■ Ar 3 (t)] , (3) 

where Afj (t) = fj (t) — fj (0) and pj (t) is the local per- 
sistence operator. That is, f>i(t) takes the value of unity 
if particle j has not moved a distance greater than some 
microscopic cutoff do, and it is zero otherwise. We use 
hats throughout this article to denote fluctuating quan- 
tities (operators). 

The usefulness of lies in the fact that the two terms 
separate mobile and immobile particles, explicitly ac- 
counting for the dynamical heterogeneity in the system. 
The expectation of each term can be simply evaluated in 
an appropriate Gaussian, or homogeneous, approxima- 
tion, leading to: 

F s (k,t) EE (F,(fc,i)} 

» P(t) + [1 - P(t)]exp(-k 2 Dt), (4) 

where P{t) is the average persistence function P(t) = 
(pj(t)) and D is the self-diffusion constant. The ap- 
proximate equality is valid [1 51 ] in the deeply supercooled 
regime of large decoupling between a- relaxation time and 
diffusion rate [20| . According to this approximation, par- 
ticle motion is a random walk with randomly distributed 
waiting times (2]| . 

Moving from two point to four point functions, we 
consider the correlator Gji(k,t) — {SFj(k,t)SFi(—k,t)) 
where 5Fj = Fj — F s was defined above. The diagonal 
(self) part is: 

Gjj{k,t) = l-F s {k,tf. (5) 



It grows monotonically with time and is always smaller 
than unity. Non-trivial spatial correlations of the dynam- 
ics appear in the off-diagonal terms j ^ I. Using Eqs. © 
and ([U, these correlations are 

G 3l (k,t)^(l-e- k2m y (6p 3 (t)6 Pl (t)} 

+ [1 - P{t)f (5 e ^-Ar 3 {t) §e -ik-Afi{t)^ (g) 

where 5pj = pj - P and <5e ife ' Af V'W = ^k-Af^t) _ e -k 2 Dt^ 
We have dropped terms that are cubic and quartic in 
the fluctuations since we expect their contributions to be 
small. 

To arrive at the four-point susceptibility, we must sum 
over j and I. The first term of ([6]) is the susceptibility of 
the persistence, and we have 

ij2( s mmt)) «iv p (t)p(i)[i-p(i)]. (?) 

Here we have assumed that excitations propa gat e 
through the system leading to movement of particles [131 ] , 
and that the average number of particles visited by a sin- 
gle excitation in time t is N p (t). This number is related 
to the dynamic correlation length of relaxed clusters of 
particles through an exponent representing their (possi- 
bly fractal) dimensionality. N p (t) is also is related to 
the average persistence: a molecule will typically persist 
after time t if there were no excitations initially in a sur- 
rounding region of mass N p (t); the probability for this 
to happen is P(t) *=a e~ cNp ^ , where c <C 1 is the average 
concentration of excitations. The calculation of the av- 
erage of the product of the persistence is analogous: for 
PjPi to be nonzero after time t a mass N p (t) has to be 
free of excitations initially around j and I. The fact that 
these two volumes may overlap gives rise to spatial cor- 
relations between pj and pi . Following similar arguments 
to those used in Ref. [1] one arrives to equation ([7]). 

While the contribution to X4(^i0 from equation (JT]) 
measures whether nearby particles relax in a correlated 
way, the second term of © measures correlations in the 
displacements of particles that have relaxed. We define 
the correlation function for particle displacements in the 
following way: 

ff( r >*) = ^ 1 {Afj(t) ■ Afi(t)) rjl=r . (8) 

where the average is conditioned on the initial separa- 
tion of the particles rji = \rj(0) — r;(0)|. This function 
g(r, t) measures the correlations between displacements 
of nearby particles. For simple Brownian motion, the 
separation of a pair of particles diffuses four times as fast 
as their centre of mass. If g{r,t) is positive, then the 
separation diffuses more slowly than this reference value; 
if it is negative then the separation diffuses faster. If 
two particles move through the system facilitated by the 
same excitation then their separation remains relatively 
small, while their centre of mass moves a long way. We 
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FIG. 2: (Color online) (a,b) Self-intermediate scattering functions F s (k,t) in the WCA mixture for wave vectors k = qo and 
k — go/4 (qo is the wave vector of the first peak in the structure factor). (c,d) Corresponding four-point susceptibilities X4(k, t). 



show the resulting positive g(r, t) for a simple kinetically 
constrained model in Fig 6 (below). 

Within the Gaussian approximation, we arrive at 



(<5e 



ifc-Ar 



-2k 2 Dt ( 2tk 2 g(r,t) 



1 



The Gaussian approximation is justified because particles 
have made many diffusive steps. The sum over particles 
at a given time t is dominated by pairs of particles whose 
initial separation coincides with the maximum of <?(r, t). 
This leads to a contribution to Xi(k, t) of 



N 



^N g {t)e- 2k " m (e 2tk ^ -1), (9) 

where g(t) = max. x [g(x,t)] obeys g(t) < D and decreases 
at large times. We have assumed that the sum over j is 
dominated by particles for which g(r,t) ~ g(t), and we 
denote the number of these particles by N g (t). The sum- 
mand is exponential in tg(r,t), so this is a good approx- 
imation when that function is large. In the excitation 
picture we only expect strong displacement correlations 
when two particles are facilitated by the same excitation, 
so we expect N g (i) < N p (t) in general. 



Putting the results of Eqs. I©-® and © together we 
arrive at the four-point susceptibility, 

X4 (fc, t) « N p (t)P(t) [1 - P(t)} (l - e- k2m 

+N g (t) [1 - P(t)f e- 2k2m [e 2tk "^) - i) 

+i- j f s (M) 2 - (io) 

Thus, correlations between particles come from persis- 
tence correlations [first term in (|10[) ] and from displace- 
ment correlations [second term in (|10p]. Both contribu- 
tions are non-monotonic in time. The persistence contri- 
bution peaks at a time i pe ak that scales as £ pC ak ~ T a , 
where r Q is the structural relaxation time. This term is 
relevant if k is large. Conversely, the second term in (fTU|) 
dominates at small k, when Dk 2 r a <C 1. The peak of this 
term occurs at £ p0 ak ~ (Dk )~ , and so increases with 
decreasing k, while its peak height decreases as <7(£ pC ak)- 

Equation (fTU|) is consistent with the data of Fig. [TJ As 
well as the overall form of Xii it contains two main pre- 
dictions. Firstly, for a given wave vector, Xi(k,t) peaks 
at a fc-dependent time £ pea k- The scaling of this time is 
the same as that of the time scale of two-point correla- 
tors [l5|. It increases with decreasing k, consistent with 
simulations. A corollary is that Xi{k,t) may display a 
non-trivial structure even for wavelengths at which the 
corresponding one-particle motion is Fickian. 
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FIG. 3: (Color online) (a) Self-intermediate scattering functions F s (k,t) in the WCA mixture at T = 0.36, for wave vectors 
k — 2qo, qo, go/2, go/4, go/8, (b) Corresponding four-point susceptibilities \4{k,t). (c) Distinct part Xi(k,t) — xT ll (k,t). (d) 
Distinct part normalized by the self term, X4,(k,t)/xT U (k,t) — 1. 



Secondly, the increase of X4(^ji) depends on the wave 
vector k: from above, we have for large k 

Xi (k,t) ~N p (t)[l-P(t)], £fc 2 i»l, (11) 

and for small k 

Xi(k,t) ~ N g (t)[l-P(t)] 2 g{t)k% DkH^l. (12) 

The two relations come from persistence (large k) and 
displacement correlations (small k). These two scaling 
predictions are different in general: we argue that this 
accounts for the variation with k at early times that was 
demonstrated in figure [TJ Within the simplest picture 
of facilitation, N p (t) is the fundamental object: it mea- 
sures the number of particles whose relaxation is facil- 
itated by a single mobility excitation. At this minimal 
level of theory, and for small times (t << t q ), we expect 
the number of particles that have relaxed to be N c N p (t), 
where N c is the number of mobility excitations; it follows 
that [1 - P(t)\ ~ N p (t). Of the N p particles facilitated 
by a single excitation, we expect a finite fraction to have 
strong displacement correlations: this finite fraction is 
Ng/Np. If this assumption holds then N g (t) and N p (t) 
will have the same scaling with time and temperature. 

Following Toninclli et al. @, we define an exponent fi 
by Xi{t) ~ & The observation of Ref. [f| is that fi ~ b, 
where b is the von-Schweidler/3-relaxation exponent of 
mode-coupling theory (MCT) [22|]. Eq. (jTTJJ) shows clearly 
that [i depends on wave vector, \x = however, MCT 



predicts a dvnamical critical point at which b is indepen- 
dent of k [22] (see nevertheless the discussion in [23l|). 
In supercooled liquids, no such phase transition is ob- 
served, and the system is always in the ergodic phase. 
Hence, all correlation lengths are finite, and the scaling 
relation /i ~ b breaks down for wave vectors smaller than 
the inverse correlation length. This effect can be investi- 
gated by direct numerical solutions of the MCT equations 
[HI HH ! but agreement with simulation is still poor for 
small wave vectors. This is consistent with the hypothesis 
that processes neglected by the mode coupling approxi- 
mation are important for the structural relaxation of su- 
percooled liquids, leading to avoidance of the dynamical 
transition; to the decoupling of diffusion and viscosity 
[25l ] ; and to non-trivial wave vector dependence of two 
and four-point correlations on large length scales. 



B. Related four-point functions 

The function 

G 4 (k,q,t) = l^(^ i (fc ) t)^(-fc.*)e < *" 1 ^ (0) -* ,(0)1 > 
ji 

generalizes x 4 (fc,i). Clearly Xi{k,t) — Gi(k,0,t). How- 
ever, G±(k,q = 0,t) is ensemble dependent (see e.g. 
(2^|). whereas lim 9 ^ Gi(k, q, t) is not. For the WCA 
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WCA, d=3, T=0.36 
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FIG. 4: (Color online) (a) Two point overlap correlators 
M a (t) in the WCA atomistic model at temperature T = 0.36. 
(b) Corresponding four point susceptibilities xm(ci, t). 



mixture and the (2)-TLG we have found that the ensem- 
ble dependence affects the absolute value of Xi(k,t) but 
not their functional trends. The differences in x&(k,t) 
between ensembles can be calculated in terms of ther- 
modynamic properties and derivatives of F s (k,t) [27| . 
Berthier et al. [26| suggest that these differences pro- 
vide reliable estimates of lim g ^o G±{k, q, t) itself. For the 
WCA mixture and (2)-TLG, in the regimes accessible to 
our simulations, we find that these differences also have 
an important k dependence, and that they only account 
for a fraction of the total value of lim ? _ ,o Gi(k, q, t). This 
is discussed in detail in the Appendix. 

An alternative (ensemble dependent) four point func- 
tion is 

XM (M) = -^^2( s Mj(a,t)6Mi(a,t)) (13) 
ji 

where SM = M — (Mi) in which Mi(a,t) is a (binary) 
operator that equals unity if \fi(t) — fj(0)| < a and zero 
otherwise. 

The operator M(a, t) was used, for example, in [1(3] ■ Its 
ensemble average, M a (t) = (Mj(a,t)), is the fraction of 
particles that have not moved beyond a distance a in time 
t, and is related to the self-intermediate scattering func- 
tion by M a (t) = (2iT)- d J M<a d d r Jd d ke^ k r F s (k,t). 

Our analysis of X4(&> t) generalizes immediately to xm- 



we write 

M 3 {a,t)=p J (t) + [l-pj(t)]e(a- |Ar 3 -(t)|), (14) 

where Q(x) is the Heaviside function. The difference 
between the operators Mj(a,t) and Fj(k,t) is that the 
phase of Fj(k,t) records the direction of motion of the 
jth particle. For small a, correlations between the direc- 
tions in which mobile particles have moved do not con- 
tribute to xm(o-i t) and we expect it to be dominated by 
the persistence correlations: 

XM (a,t) « f(Dt/a 2 )N p (t)P(t)[l-P(t)} (15) 

where f(x) = (4irDt)- d [f H<a d d r e -k| 2 /(£>*)]2. Thc 
structure is the same as the first term of (fTU|) with 
f(Dt/a 2 ) playing the part of (1 - er Dk2t ) 2 . We note 
also that XM(a,t) approaches the contribution of (0 as 
a gets small, as expected. We compare this prediction 
with atomistic simulations in section [TTTl see below. 



III. ATOMISTIC SIMULATIONS 

For a continuous atomistic representation of a super- 
cooled liquid, we carried out extensive molecular dynam- 
ics simulations of a symmetric WCA mixture [16| . It 
is a mixture of two particle species A and B in three 
spatial dimensions. The potential energy is the sum of 
the pairwise interactions between two particles of species 
a and /3, V a p(r) = 4e (a Q p/r) 12 - {<r al3 /r) 6 + 1/4 if 

r < 2 1 / 6 er Q/ 3, and V a p(r) = otherwise. This is the 
reference potential of the WCA theory [17] consisting 
of the repulsive part of the Lennard- Jones interaction. 
Following [nj we choose: ctaa — I, &bb — 5/6, 
°ab = {p aa + gbb)I% m B = m A /2 = I, and s = 1. 
Lengths, times and temperatures are given in units of 
ctaa, \J fnBC AA /e, and e/ks, respectively; we use qo to 
denote the wave vector of the first peak in the struc- 
ture factor. In our simulations the total number of par- 
ticles was N = 8000, with N A = N B = 4000. The use 
of the WCA reference potential makes this system com- 
putationally more efficient to simulate than the original 
Lennard- Jones one of Ref. [28[ . A detailed study of the 
dynamics in the supercooled reg ime of this WCA mixture 
will be presented elsewhere [29(. The molecular dynam- 
ics simulations of this section conserve energy, and we 
use a range of energies at each temperature. This allows 
estimates of both canonical and microcanonical suscep- 
tibilities. Details are given in the Appendix; the data of 
Figs. [2]|4]is for the canonical susceptibility. 

Figure [2] shows the self-intermediate scattering func- 
tion F s (k,t) at various temperatures T — 0.5 to 0.36, 
for the wave vector of the first peak of the structure fac- 
tor k = go, and for a smaller wavevector k = go/4. At 
the lowest temperature shown, T — 0.36, the system is 
clearly in the supercooled regime, dynamics is heteroge- 
neous, and the self-diffusion constant exceeds the value 



expected by the Stokes-Einstein relation by over an order 
of magnitude [29(. The four-point susceptibilities have 
the expected behaviour, becoming larger and peaking 
later the lower the temperature, as shown in the lower 
panels of figure [2] 

Figure [3] concentrates on the lowest temperature we 
simulated, T = 0.36. It shows the two point function 
F s (k,t) for various wave vectors k. The growth of the 
distinct part of xa towards its peak can be fitted by a 
power of t. The fitted exponent is smaller than 1 for 
wave vectors near go, but closer to 1.5 for small k. The 
former result coincides with what was observed in similar 
systems H, but the latter was not anticipated before. 
Interestingly, a similar k dependence has been observed 
experimentally in a sheared granular material 

In section |nl we assumed that XM(<2,i) is dominated 
by persistence correlations for small a, where displace- 
ment correlations are unimportant. At small times, the 
effect of increasing a is to reduce the susceptibility. This 
is analogous to the effect of decreasing k in Xi(k, t)- On 
the other hand, the k dependence of the peak in Xi(k, t) 
arises from the displacement correlations and will be ab- 
sent (or at least much weaker) in XM(o,,t). Both these 
predictions are consistent with the data of Fig. [4] cal- 
culated in the atomistic system of section ITO1 However, 
we note that the relative range of a over which we have 
measured Xm («, t) is smaller than the range of k used 
for X4(k,t). It may be that particle displacement corre- 
lations will be important at larger a. This effect can be 
estimated within the framework of section [II] It appears 
as an extra term in (11511. 
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IV. KINETIC ALLY CONSTRAINED MODELS 



FIG. 5: (Color online) (a) Self-intermediate scattering func- 
tion F 3 (k,t) of particles in the (2)-TLG model, at density 
p — 0.8, for wave vectors k — n, tt/2, tt/8, n/16, tt/32, 7r/64. 
(b) Susceptibilities x±{k,t) for the self-correlators of the top 
panel. For large k, Xi grows with an exponent close to 0.7, 
for small k with an exponent close to 1.5. The inset shows 
the normalized function X4/xl c ' f — 1 in the same scale as the 
main panel. 



A. Kinetically constrained lattice gas 

The kinetically constrained lattice gases of Jackie and 
Kronig [l8[ provide simple caricatures of supercooled liq- 
uids. We consider the (2)-TLG in which hard core parti- 
cles move on a triangular lattice: movement between sites 
and i and j is allowed only if both the mutual neighbours 
of sites i and j are empty. Relaxation in the model in- 
volves strongly co-operative motion and the relaxation 
time increases very rapidly with increasing density. 

We show correlation functions in Fig. [5] We plot the 
self-intermediate scattering functions for particles at a 
high density p — 0.8; we use wave vectors from k = tt to 
k = 7r/64 (for details see Ref. [3(|). We note that the ex- 
ponents governing the growth of X4(k,t) in the (2)-TLG 
are rather similar to those in the atomistic system. These 
four-point functions are measured at constant density. 

B. Fredrickson-Andersen (FA) model 

The one-spin facilitated Fredrickson-Andersen (FA) 
model represents the extreme of coarse-graining in which 



one conjectures that the only variables relevant for het- 
erogeneity are binary labels for regions in which mobil- 
ity is present. It is the simplest kinetically constrained 
model displaying both dynamic heterogeneity [J] and de- 
coupling [13T ], The system is a one dimensional chain 
in which mobile regions are represented by up spins and 
immobile regions by down spins. The energy of the sys- 
tem is simply the number of up spins: spins flip with 
Metropolis rates if and only if at least one neighbouring 
spin is up. 

In order to consider particle motion in such a sys- 
tem, we couple probe particles to the spins as in [13j. 
Probes can hop only between adjacent mobile sites and 
such moves are attempted with rate unity. Probes do not 
interact with one another (more than one probe may oc- 
cupy a single site). The stationary distribution of probes 
is uniform and uncorrelated with the spin variables. 

Figure [5] (a) shows the self-intermediate scattering 
functions of probe molecules in the FA model at a low 
temperature T — 0.2 for various wave vectors (the ensem- 
ble is again that in which the probe density is constant). 
For large wave vector the self-correlator is stretched and 
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FIG. 6: (Color online) (a) Self-intermediate scattering 
function F s (k,t) of probes embedded in an FA model in 
d = 1, at temperature T = 0.5, for wave vectors k — 
7r, 7r/ll, 7r/31, 7r/51, 7r/71, 7r/91. The inset shows the suscep- 
tibility of the persistence given in Eq. (0 (dot-dashed) and 
the quantity N p (t) (full), which grows approximately as t 1/l2 
(dashed), (b) Susceptibilities x^(k,t) for the self-correlators 
of the top panel. We show the non-trivial part, with the self- 
term X4° lf {k, t) removed. Xi f° r k =ir follows the contribution 
of Eq. (J7J) (dot-dashed). For large k, Xi grows with an expo- 
nent close to 1, for small k with an exponent close to 1.5. The 
inset shows the function g(r,t), see Eq. (JS|, for inter particle 
distances r = 0, 1, 2, in the same temporal scale as the main 
panel. 



tracks the persistence function, while for small enough 
wave vectors self-correlators become exponential, accord- 
ing to Eq. ([3]) (see Ref. [l5[ for details). The inset shows 
the time dependence of N p (t), as extracted from the 
susceptibility of the persistence, see Eq. ([7]). Figure [6] 
(b) shows the corresponding susceptibilities X4(k,t). We 
note that the exponent governing the increase of X4(/c, t) 
is near to unity for large k, as reported in [8J; however 
the exponent at small k is closer to 1.5. In the notation 
of mi) we have N p (t) ~ [1-P(t)] ~ t 1 / 2 at large k which 
accounts for the exponent of unity in that regime. The 
inset to the bottom panel shows the function g(r,t), for 
r = 0, 1, 2. As t increases the distance r over which these 
correlations are significant also increases. This is consis- 
tent with our assumption of an increasing function N g (t) 
in section [TT1 
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APPENDIX: ENSEMBLE DEPENDENCE OF X 4 

Recently, Berthier et al. !26j suggested that the dif- 
ference between Xi(k,t) in two different ensembles pro- 
vides a good estimate of the value of Xi{k,t) itself. In 
this Appendix we discuss the ensemble dependence of the 
four point susceptibility in the models that we consider. 
Not surprisingly, we find that the difference in Xi{k,t) 
between ensembles has a significant wavelength depen- 
dence. At the lowest temperature we simulated for the 
WCA mixture, and the highest density for the (2)-TLG, 
we find that the method of [2(| does not give accurate 
numerical estimates of the susceptibility, although it does 
give the correct order of magnitude for Xi{k,t) at times 
for which that function is near its maximum, the estimate 
working better for smaller k than for larger k. It would 
be interesting to see whether the estimate of [26( becomes 
better or worse with decreasing temperature, but this is 
beyond the scope of this appendix. 

In terms of equation ([1]), the ensemble dependence of 
Xi(k, t) is contained in the definition of the average re- 
quired by the angle brackets. This average is over any 
stochasticity in the dynamics, and over an initial condi- 
tion that is specified by the choice of ensemble. The key 
point is that X4,(k,t) measures fluctuations in an observ- 
able that couples with equal weight to every particle in 
the system. So the presence of a global constraint (on 
the total energy of the system, for example) suppresses 
these fluctuations. By contrast, the quantity G±(q,k,t) 
measures fluctuations on finite length scales, as long as 
k and q are both finite. In that case, taking the ther- 
modynamic limit leads to a value of G±(q, k,t) that is 
independent of global constraints. 

The effect of global constraints on Xi{k, t) can be mea- 
sured directly in the (2)-TLG, by considering the suscep- 
tibility in both canonical and grand canonical ensembles. 

We denote the canonical susceptibility by X^i^tt), 
which is defined as in (P), with the averages over the 
dynamics, and over all initial conditions with exactly N 
particles. For grand canonical systems, we promote N to 
a fluctuating quantity (operator), and write: 

N 2 

xi fl) (k,t) = (N)- 1 (\^2(F j (k,t)-F s (k,t))\ ). (A.l) 
where the average is over the dynamics, and over initial 
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FIG. 7: (a) Figure showing Xi(k, t) in the (2)-TLG at p = 0.77 
and k — tv/2 in both canonical and grand canonical en- 
sembles, together with the difference term Axi '{k,t) = 
p 2 (l - p)[dF s {k,t)/dp] 2 . The contribution of Axi M> (M) to 
X4*'(fc, t) is comparable to that of x^\k,t). (Inset, a) Il- 
lustration that the grand canonical susceptibility satisfies the 
relation (|A.2[I . The scales are those of the main panel, (b) 
Similar data in the (2)-TLG at k = For smaller wave 

vectors, such as these, the difference term dominates the peak 
of the grand canonical susceptibility. (Inset, b) Two contri- 
butions to the grand canonical susceptibility: see (|A.4|) . At 
early times, the susceptibility is dominated by the ensemble 
independent sine part, while the peak is dominated by the 
cosine part. The scales are again those of the main panel. 



conditions with all particle numbers, with weights set by 
the chemical potential pi. 

These susceptibilities are related to each other by [27j 



xi M) (M) 
Axi M) (M) 



xi N \k,t)+A X r(k,t), 

V Pit? II, +m 2 

P V{(6pf) 



.00 
dF s (k,t) 



dp 



(A.2) 
(A.3) 



where p is the density and V the volume of the sys- 
tem; the average is grand canonical. The density fluc- 
tuations of the (2)-TLG are those of a lattice gas, so 
that V({Sp) 2 ) = p(l - p). 

The relative sizes of the two terms on the right hand 
side of (|A.2|) depend on the strength of the coupling of the 
dynamics to density fluctuations in the system. Berthicr 
et al. argued [26[ that the second term tends to domi- 
nate the grand canonical susceptibility in athermal glassy 
systems. Since the thermodynamic properties of the (2)- 



FIG. 8: (a) Figure showing four point susceptibilities in the 
WCA mixture at T = 0.36 and k = qo, and the difference term 
Ax4 T ^(fc,t). The difference term is significantly smaller than 
the canonical susceptibility, (b) Similar data at k — qo/8. 
In this case the difference term exceeds the microcanonical 
susceptibility. The data is consistent with the exact rela- 
tion (|A-10J> - 



TLG system are trivial, Xi{k, t) can be measured in both 
ensembles, and the prediction for the difference can also 
be evaluated. Results are shown in figure at a high 
density p = 0.77. The difference between ensembles is 
consistent with the prediction of (|A.2[) . We find that the 

two contributions to t) are °f similar size for large 

wave vectors, while the second term dominates at small 
wave vectors. In both cases, the form of the second term 
is similar to that of x£ (k, t) for times near the peak of 
the susceptibility. 

Our definition of Xi{k,t) includes both the sine and 
cosine parts of the operator Fj(k,t) — cos[fc • Afj(t)] + 
isin[fe • Afj(i)], leading to two contributions to X4 
That is, 



X4 (M) = X4 os (M) + X4 in (M), (A.4) 



with 



xT(k,t) = (N)-^{j2co S [k- AfM -F B (k,t)} 

(A.5) 

Xf n (k,t) = (^^({^sinlfc-Ar^t)]} 2 ), (A.6) 
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The fluctuations in the sine part are ensemble indepen- 
dent, since the average of this part is zero. Hence, the 
difference term Ax± (k,t) is bounded above by the co- 
sine part of (k, t); it follows that the difference term is 
qualitatively different from the total susceptibility when 
the total susceptibility is dominated by the sine part. 
This is the case at small wave vectors and small times, 
as shown in the inset of the lower panel of figure [7] 

We can make a similar analysis of the WCA mixture. 
For each temperature, the data of section IIIII was ob- 
tained using microcanonical simulations at eight differ- 
ent values of the energy, sampled from the appropriate 
canonical ensemble. This allows estimates of both mi- 
crocanonical and canonical susceptibilities, which we de- 

( E) (T) 

note by X4 (k,t) and \a \k, t) respectively. The micro- 
canonical susceptibility for a given energy is given by |T]), 
with microcanonical averages throughout. In particular, 
we must use 

microcanonical; 

(A.7) 

where the average involves data at a single energy. Hav- 

( E) 

ing used this prescription to calculate x\ (k, t) for each 

( E) 

energy, we then average xk {k, v over the eight repre- 
sentative energies at each temperature. Since the fluctu- 

( E\ 

ations in x\ (k,t) itself are small, the resulting average 
is an unbiased estimate of the microcanonical suscepti- 
bility at the (canonical) average energy associated with 
that temperature. For the canonical susceptibility, we 
use ([1]), with canonical averages throughout, including 

canonical • 

(A.8) 

We estimate canonical averages by averaging over the 
data from the eight energies that are representative of 
the relevant temperature. 

The canonical and microcanonical susceptibilities sat- 
isfy [13 

xi T) (M) = xi B) (M) + Axi T) (M) (A.9) 

A x i T) (fc,<) = (k B T 2 /c v )(dF s (k,t)/dT) 2 (A.10) 



where cy is the specific heat per particle at constant vol- 
ume. The behavior of cy that we have computed for the 
WCA mixture is unremarkable. It rises slowly with low- 
ering temperature, indicative of mean-square potential 
energy and kinetic energy fluctuations that are virtually 
independent of temperature over the lowest temperatures 
we have considered. 

We show results for both susceptibilities in figure [5J 
together with an estimate of AX4 (k, t). [In evaluating 
^■x\ {k, t), the derivative and specific heat are evaluated 
by a finite difference analysis of the temperature depen- 
dence of the the mean energy and F s (k,t)j. The data is 
at the low temperature T = 0.36. At large k, x\ (k,t) 
is not saturated by the difference term Ax\ (k,t). In- 
stead the microcanonical contribution is larger. However, 
at smaller k, the difference term exceeds the microcanon- 
ical part for times near the peak. Just as in the (2)-TLG, 
the difference term becomes more significant at smaller 
k, and at later times. Since the size of the difference 
term reflects coupling between energy fluctuations and 
dynamics, it would appear that motion over large dis- 
tances couples strongly to the local value of the energy 
while local rearrangements do not. 

Very recently, Szamel and Flenner (3lj mea- 
sured the contribution of global energy fluctua- 
tions to the canonical susceptibility Ax 4 (k,t) = 
(k B T 2 /c v ° t )(dF s (k,t)/dT) 2 , in a model glass-former 
with Brownian dynamics (here, Cy 0t is the potential en- 
ergy contribution to the heat capacity, at constant vol- 
ume). Since Brownian dynamics are non-conservative, 

Ax\ [k, t) does not represent a bound on the canoni- 
cal X4(k, t). However, the results of [3l| show that while 
^Xi [k, t) and the unconstrained susceptibility are not 
equal, they are of the same order of magnitude for large 
k. We find a similar situation in the (2)-TLG; in the 

WCA mixture then the difference term Ax\ (k,t) is 
much smaller than the canonical susceptibility at these 
wave vectors. 
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